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^ ' We present an exploratory QCD analysis of the neutrino structure functions in 

OO ■ charged current DIS using the color dipole formalism. The dipole cross sections are 

taken from recent phenomenological/theoretical studies in deep inelastic inclusive 
Q^ I production. The theoretical predictions are compared to the available experimental 

■ results in the small-x region, which has never been considered so far. 
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1 Introduction 



X 

a I The interaction of high energy neutrinos on hadron targets are an outstanding 

probe to test Quantum Chromodynamics (QCD) and understanding the par- 
ton properties of hadron structure. The several combinations of neutrino and 
anti-neutrino scattering data can be used to determine the structure functions, 
which constrain the valence, sea and gluon parton distributions in the nucle- 
ons/nuclei. The neutrino structure functions are needed for computing the to- 
tal neutrino-hadron cross section, which plays an important role in high energy 
cosmic rays studies and in astroparticle physics [1]. The differential cross sec- 
tion for the neutrino-nucleon charged current process vi {ui)-\-N —>■ l~ (/"*") +X, 
in terms of the Lorentz invariant structure functions -^2''^? 2xF^'^ and xF^^ 
are [2], 



dxdy 7T 
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where Gp is the weak Fermi couphng constant, mjv is the nucleon mass. E,y is 
the incident neutrino energy. is the square of the four-momentum transfer 
to the nucleon. The variable y = Ehad/Ei, is the fractional energy transferred 
to the hadronic vertex with Ehad being the measured hadronic energy, and 
X — /2mNEi,y is the Bjorken scaling variable (fractional momentum carried 
by the struck quark). 



Similarly to the charged- lepton DIS, the deep inelastic neutrino scattering 
is also used to investigate the structure of nucleons and nuclei. In the lead- 
ing order quark-parton model (the QCD collinear approach), the structure 
function F2 is the singlet distribution, F^^ oc xq^ = xJ2{<1 + q): the sum 
of momentum densities of all interacting quarks constituents, and xF^ is the 
non-singlet distribution, xF^^ oc xq^^ = x^{q — q) = xuv + xdy, the valence 
quark momentum density. These relations are further modified by higher-order 
QCD corrections. Currently, the theoretical description of experimental data 
on neutrino DIS is reasonable (for a very recent investigation, see Ref. [3]). 
The main theory uncertainties are the role played by nuclear shadowing in 
contrast with lepton- charged DIS and a correct understanding of the low 
limit. The first uncertainty can be better addressed with the future precise 
data from MlNERz/A [4]. However, nuclear effects are taken into account by 
using the nuclear ratios R = F2/AF2 extracted from lepton- nucleus DIS, 
which it could be different for the neutrino-nucleus case. The low-Q^ region 
can not be addressed within the pQCD quark-parton model as a hard mo- 
mentum scale Qq>1 — 2 GeV^ is required in order to perform perturbative 
expansion. 



In this work we present a determination of the small-x structure functions for 
neutrino- nucleus within the color dipole formalism [5]. This approach allows 
for a simple implementation of shadowing corrections [6] in neutrino-nuclei in- 
teractions. This paper is organized as follows. In Sec. 2, the structure function 
F2^ is investigated within the color dipole picture at small-a: region, em- 
ploying recent phcnomenological parton saturation models. It is shown that 
small-a; data exhibit geometric scaling property, which has important conse- 
quences for ultra- high energy neutrino phenomenology. In Sec. 3, the structure 
function xF^ and the quantity AxFs are addressed. The latter one provides 
a determination of the strange-sea parton distribution through charm pro- 
duction in charged-currcnt neutrino DIS. Finally, we also analyze the nuclear 
ratios R2 = F^^/AF^^ and R2 = xF^'^/AxF^'^ . In the last section we present 
comments and conclusions. 
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2 Neutrino Structure Function F^^{x, Q^) 



We focus here in the high energy regime, which one translates into small-x 
kincmatical region. At this domain a quite successful framework to describe 
QCD interactions is provided by the color dipole formalism [5], which allows an 
all twist computation (in contrast with the usual leading twist approximation) 
of the structure functions. The physical picture of the interaction is the deep 
inelastic scattering viewed as the result of the interaction of a color qq color 
dipole, in which the gauge boson fluctuates into, with the nucleon target. The 
interaction is modeled via the dipole-target cross section, whereas the boson 
fluctuation in a color dipole is given by the corresponding wave function. The 
charged current (CC) DIS structure functions [7-9] are related to the cross 
section for scattering of transversely and longitudinally polarized bosons. 
That is, 

= ^ / £ \'^T,t {z, r, Q') \' a^, (x, r) , (1) 

where r denotes the transverse size of the color dipole, z the longitudinal mo- 
mentum fraction carried by a quark and ■j/'t^l are the light-cone wavefunctions 
for (virtual) charged gauge bosons with transverse or longitudinal polariza- 
tions. The small- a; neutrino structure function F!^^ is computed from expres- 
sions above taking F2 = Ft + Fl. Explicit expressions for the wave functions 
squared can be found at Refs. [7-9]. In what follows we consider four quark 
flavors {u, d, s, c) with masses m/. The color dipoles contributing to Cabibbo 
favored transitions are ud{du), cs{sc) for CC interactions. 

The dipole hadron cross section aaip contains all information about the target 
and the strong interaction physics. In the present study, we consider analytical 
expressions for the dipole cross section, with particular interest for those ones 
presenting scaling behavior. Namely, one has adip oc {r'^QsatV dipole sizes 
~ 1/Osat ^'^d where (1 — 7) is the effective anomalous dimension. The so- 
called saturation scale Qsat oc x^^'^ deflnes the onset of the parton saturation 
effects. In what follows one takes the phenomenological parameterizations: 
(a) Golec-Biernat-Wiisthoff model (GBW) [10] and (b) Itakura-Iancu-Munier 
(IIM) model [11]. Both models are able to describe experimental data on in- 
clusive and diffractivc deep inelastic ep scattering at small-x. We quote the 
original papers for details on the parameterizations and determination of their 
phenomenological parameters. We call attention to Ref. [12], where the im- 
plications of saturation in ultrahigh energy neutrino cross section has been 
flrst discussed within the color dipole picture. We will show that the results 
are not strongly sensitive to a particular choice of model. Here, we use an 
effective light quark mass, rrif = 0.14 GeV and the charm mass is set to be 
rric — 1.5 GeV. In addition, as the color dipole models are suitable in the 
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Fig. 1. The scaling behavior of crwNixjQ"^) as a function of scahng variable r (see 
text). 

region below x = 0.01 and the large x limit still needs a consistent treatment, 
we supplement the dipole cross sections with a threshold factor (1 — 3;)"'"^'^'== 
('^thres = 5 (7) for uumbcr of flavors nj = 3 (4)). 

The extension of the approach to consider nuclei targets we take the Glauber- 
Gribov picture [6], without any new parameter. In this approach, the nuclear 
version is obtained replacing the dipole-nucleon cross section by the nuclear 
one. 
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where b is the impact parameter of the center of the dipole relative to the 
center of the nucleus and the integrand gives the total dipole-nucleus cross 
section for a fixed impact parameter. The nuclear profile function is labeled 
hjTAib) [13]. 

A basic property of the saturation physics is the geometric scaling. It means 
that the total 7*p cross section at large energies is not a function of the two 
independent variables x and Q, but is rather a function of the single variable 
Tp = Q"^ /Ql^^{x) as shown in Ref. [14]. That is, a^*p{x, Q"^) = a-f*p{Tp). In Refs. 
[15, 16] it was shown that the geometric scaling observed in experimental data 
can be understood theoretically in the context of non-linear QCD evolution 
with fixed and running coupling. Recently, the high energy /^p, pA and AA 
collisions have been related through geometric scaling [17]. Within the color 
dipole picture and making use of a rescaling of the impact parameter of the 7*/i 
cross section in terms of hadronic target radius Rh, the nuclear dependence 
of the 7*^4 cross section is absorbed in the A-dependence of the saturation 
scale via geometric scaling. The relation reads as (Jtot^iTA) = t^A(^tot {'^a)-, 
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Fig. 2. The structure function F!^^ {x,Q'^) as a function of boson virtuality. The 
long-dashed curve corresponds to the geometric scahng result for the neutrino-nu- 
cleus interactions, the dot-dashed curve to the valence content and the solid one 
is the sum of both contributions. Experimental data from CCFR Collaboration 
[19,20]. 

where ka = {Ra/RpY- The nuclear saturation scale was assumed to rise with 
the quotient of the transverse parton densities to the power A and Ra is the 
nuclear radius, Qsat.A — (^/'^a)^ QLt.p- The functional shape for the pho- 
toabsortion cross section, a^*p (Tp), has been considered based on theoretical 
studies [17]. 

The geometric scaling property has direct consequences on the computation 
of small- a: neutrino structure functions. It has been shown in Ref. [18] that the 
charged (CC) and neutral (NC) current structure functions are described by 
the same mathematical expressions as the proton structure function up to a 
different coupling of the electroweak bosons. Thus, using such a property one 
has for the neutrino-nuclei collisions [18], 



'tot 



X. 



Q'; A) = ao 



'A, 



(3) 



where is the QED constant coupling, e/ is the sum of electric charge of 
the quarks of flavor /. The quantity •je is the Euler constant and T (0, /?) the 
incomplete Gamma function. The parameters a and b were obtained from a 
fit to the small-x ep DESY-HERA data, with the overall normalization fixed 
by ao = 40.56 /xb [17]. 
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It is important to investigate whether this geometric scaUng pattern is exhib- 
ited by the smaU-x neutrino deep inelastic data. In order to do so, we take 
the datasets for the structure function F2 with the kinematical cut x < 0.035 
and all Q'^ [19-21]. These experiments measure the differential cross sections 
for deep inelastic i^^-Fe and P^^-Fe scattering and the structure functions are 
then extracted. In Fig. 1, we plot the quantity awN — ^F^^ as a function 
of the scaling variable r = Q"^ /Ql^^{x). The saturation scale squared reads as 
Qsat(^) ~ (3 ■ 10^'^/a;)°-^^^ [10]. The geometric scaling behavior is present with 
a smooth spread around the scaling curve. The scaling curve (solid line) is 
obtained using Eq. (3). This result is useful and can be investigated in more 
detail using precise measurements in future experiments. In the present kine- 
matical window, the color dipole formalism (and geometric scaling property) 
is in the limit of its validity. However, the present result shows that it gives a 
reasonable phenomenological description of the limit case. 

Now, we will compare the color dipole prediction against the structure function 
= j:^(JwN- This is presented in Fig. 2. We use the experimental datasets 
of the CCFR Collaboration [19, 20], where filled circles correspond to points in 
Ref. [20] and triangles up correspond to points in Ref. [19]. The long-dashed 
curve is obtained using scaling expression Eq. 3) for the boson-nucleus cross 
section. Nuclear effects are taken into account through the nuclear saturation 
scale. The calculation produces a suitable description at small-x, despite the 
data points lying at the expected validity region of the color dipole approach, 
X < 10^^. For completeness, we have added the valence content to F2. To this 
aim, the following parameterisation has been considered [22], 



which it is a Regge inspired model. Above, n(Q'^) = § (l + ^^)- The pa- 
rameters are taken from Ref. [22], with B,y = 2.695, fi, = 0.595, aR — 0.425, 
c = 3.5489 GeV^ and b = 0.6452 GeV^. The valence contribution is repre- 
sented by the dot-dashed curve, which is subleading at small- a;, but improves 
the overall description. The total contribution (valence -|- sea color dipole) is 
given by the solid curve. It is worth to mention that we have not tunned the 
original parameters of the scaling model and of the valence parametrization. 

A short comment is in order here. We found that the complete numerical cal- 
culation for the structure function F2, using Eqs. (1) and nuclear corrections 
from Eq. (2) produces results compatible with the geometric scaling calcu- 
lation above within a few percents. This fact corroborates the advantage of 
using the fast scaling parametrization, which can be promptly computed for 
any nuclei. This is possible because the parameters of the scaling curve suit- 
ably absorb the details of the charm contribution, which turn out the complete 
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Fig. 3. The structure function xF^^ (xjQ"^) as a function of boson virtuality (see 
text). 

calculation somewhat more involved. In the next sections, we will use the full 
color dipole calculation as the structure function xF^ is strongly sensitive to 
the charm quark contribution. 



3 The quantities xF^'^ , AxF^'^ and the nuclear ratios 



Lets now compute the structure functions xF^ within the color dipole formal- 
ism. We concentrate on the interaction of the cs color dipole of size r with the 
target hadron which is described by the beam- and flavor-independent color 
dipole cross section a^ip- In the infinitum momentum frame, this is equivalent 
to the M/^^-gluon fusion process, + 9 cs (cs). The analysis for charged 
current DIS has been addressed in Refs. [23,24], where the left-right asym- 
metry of diffractive interactions of electroweak bosons of different helicity is 
discussed. There, the relevant light-cone wavefunctions have been evaluated. 
The contribution of excitation of open charm/strangeness to the hadron ab- 
sorption cross section for left-handed (L) and right-handed (i?) VT-boson of 
virtuality Q^, is given by [5], 



<Jl,r{x,Q') = d'r dzYl \^l]f{z,r,Q')\'a,,p{x,r), (5) 

° Ai,A2 

where "^^'^^{z, r, Q"^) is the light-cone wavefunction of the cs state with the c 
quark carrying fraction z of the light-cone momentum and s with momen- 
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turn fraction 1—z. The c- and s-quark hclicitics arc Ai ~ ±1/2 and A2 = il/2, 
respectively. The diagonal elements of density matrix are given by, 



AN, 



Ai,A2 



(27r)2 
AN, 



c 2 
Z 



m\K'^{er)+e'Ki{er) 



Y: K"'- (^^^^'f = (1 - zf [ml Kl{er) + Kl{er) 



(6) 
(7) 



In the expressions given by Eqs. (6)-(7), one uses the notation — z{l — 
z)Q^ + (1 — z)m1 + -2m?, where the quark and antiquark masses are iriq and 
TJiq, respectively. The corresponding expressions for W~ boson are obtained 
by replacing niq ^ rriq. It should be noticed the strong left-right asymmetry 
referred above. 

The structure function of deep inelastic neutrino-nucleon xF^ can be defined 
in terms of and ai of Eq. (5) in the following usual way. 



xFf (x, Q') = £j [a^x, Q') - an{x, Q')] ■ (8) 

where the expression can be interpreted in terms of parton densities as being 
the sea-quark component of 2:^3. It corresponds to the excitation of the cs 
state in the process W^g — > cs, with xF^ differing from zero due to the strong 
left- right asymmetry of the light-cone \cs) Fock state. For values of Bjorken 
variable not so small, xF^ contains important valence quark contribution. The 
valence term, xq^aU is the same for both vN and vN structure functions of an 
iso-scalar nucleon. The sea-quark {xqsea) term in the xF^^ has opposite sign 
for xF^^ , leading to xF^'-"^^ = xq^ai ± xqsea- 

A direct comparison of result in Eq. (8) with experimental data is somewhat 
difficult. The reason is that structure functions obtained from neutrino scat- 
tering experiments are usually extracted from the sum and the difference of 
the neutrino and anti-neutrino y-dependent differential cross sections, respec- 
tively. That is, xF^ is determined by the average \{xF^^' + xF^^). In Ref. [23], 
the authors intend to compare the result of Eq. (8) to CCFR data [19]. How- 
ever, this procedure is questionable since the cs(sc) component disappears in 
the sum F^ + F^ = 2xqyai- This comparison is shown in Fig. 3, where the 
structure function xF^^{x, Q"^) is shown as a function of Q'^. The theoretical 
curves correspond to the results using the GBW (dashed-line) and the IIM 
(dot-dashed line) models for dipole cross section, respectively. The total con- 
tribution (solid) includes the valence quark contribution through the above 
procedure. The calculation takes into account nuclear shadowing for the iron 
nucleus via Glauber-Gribov formalism, where the dipole-nucleus cross section 
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Fig. 4. The structure function AxF^^ as a function of boson virtuality (see text). 

is computed using Eq. (2). The distinct models for the dipole cross section 
give similar results for the range of virtuality considered in the plots. We have 
checked our results with those in Ref. [23] and it was found they are similar 
despite the different dipole cross section used. The main reason is that current 
experimental results are probing mostly the color transparency domain in the 
dipole cross sections. Future measurements at smaller x and heavy nuclei, as 
proposed in the Minerva experiment [4], will allow to disentangle the models 
concerning the parton saturation effects. 

The neutrino-antineutrino difference xF^ — xF^ provides a determination of 
the sea (strange) density. In the parton model, one has xFg^^ = xq^ai — 
2xc{x) + 2xs{x) and xF^^ = xq^ai + 2xc{x) — 2xs{x). Therefore, the neutrino- 
antineutrino difference effectively measures the strange density, since the charm 
contribution is small in the kinematical region measured by current experi- 
ments. Assuming s{x) = s{x) and c{x) = c{x) one obtains, 

AxFs = xF^^ - xFl^ = 2xqsea = 4x [s{x) - c{x)] . (9) 

Here, some comments are in order. Taking into account that our calculation 
in Eq. (8) corresponds to sea-quark content oi xF^, then AxF^, = 2xF^^ . Our 
calculation is equivalent to the cs component of the structure function given 
by the M/^-gluon fusion term at order a^, which reads as, 

Ff{W^g cs) = (^) £ ^g{z, /.^) C, Q^) , (10) 
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where a — 1 + (m^ + ni'D/Q'^ and the Wilson coefficient C3 represents the 
W^g cs cross section. The total cs contribution to F3 is the sum of the 
quark excitation term, taken at the factorization scale fi^ = rn^ (near thresh- 
old), and the gluon-fusion term above. Namely, F^^ (x, Q'^) = 2 [s{xc, /x^) — c{x, + 
F^^{W'^g — > cs), with the slow-rescaling variable Xc — x[l + {rric/Q'^)] and 
similar expression for F^^ (x, Q'^). 

In Fig. 4 the quantity AxF^ as a function of at fixed x is shown in com- 
parison with the CCFR result obtained from u^Fe and P^Fe differential cross 
section [25]. The theoretical curve is obtained from Eq. (9) using the IIM 
dipole cross section and Glauber-Gribov shadowing corrections. The agree- 
ment is good and the quark excitation contribution, described above, has not 
been added. This additional piece should improve the description. 

As a final analysis, we consider the nuclear ratios for the neutrino structure 
functions. We are confident in the reliability of these calculations using the 
color dipole approach as they were successfully tested against charged lepton 
scattering data [6]. First, we compute the ratio R2{x,Q'^) = F2'^ / {AF!^'^). 
This is shown in Fig. 5 function of for fixed x (IQ-^ <x< lO'^). 

For our purpose, an iron nuclei {A = 56) is considered and valence content 
is disregarded. A strong dependence is observed. For intermediate virtu- 
alities, the ratio ranges on 0.75 < R2 < 0.85 and decreases with lower x. In 
the current region of neutrino data, the nuclear correction is of order 15-20 
%. This dependence on virtuality is different of the observed flat behavior 
present in calculation using the DGLAP collinear approach [3]. In Fig. 5-b, it 
is shown the nuclear ratio Rz{x, Q^) = xFJ^^/^AxFJ^'^) as a function at at 
fixed X. This corresponds to nuclear shadowing in sea content of xF^. A weak 
dependence on Q"^ is observed, whereas it decreases at lower x. An interesting 
point it is the different shadowing contribution for xF^ in contrast with F2. 
This feature has already been observed in the comparison of collinear QCD 
approach with neutrino data [3] and corroborates the dynamical shadowing 
correction considered here. 



4 Comments and Conclusions 

As a summary, an analysis of small-x neutrino-nucleus DIS is performed within 
the color dipole formalism. The structure functions F2^ , xF^^ and the quan- 
tity AxFg ^ are calculated and compared with the experimental data from 
CCFR and NuTeV by employing phenomenological parameterizations for the 
dipole cross section which successfully describe small- a; inclusive and diffrac- 
tive ep DIS data. Nuclear shadowing is taking into account through Glauber- 
Gribov formalism. It is found that small-a; data show geometric scaling prop- 
erty for the boson-hadron cross section as a function of the scaling variable r. 
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Fig. 5. The nuclear ratios (a) i?2 = F^I{AF2) and (b) i?3 
function of virtuality at fixed x (see text). 



(b) 

as a 



The structure function F2 is in agreement with the phenomenological imple- 
mentation using the saturation models at the small-x region. The structure 
functions xF^^ is also discussed in detail. The sea content, described by the 
quantity AxF^^, is well described and the addition of quark excitation term 
should improve it. Although the results presented here are compelling, further 
investigations are requested. In particular, measurements of neutrino-nucleus 
structure function in smaller values of x than the currently measured in the 
accelerator experiments. The present results also confirms the robustness of 
the color dipole formalism to describe the total neutrino cross section of ultra- 
high energy neutrinos. Finally, we predict the nuclear ratios R2 and and 
single out the size of nuclear effects in each case. 
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